"school 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

EFFICIENT  COMPUTATION  OF  THE  K-T 
RELIABILITY  OF  DIRECTED  ACYCLIC 

ERMINAL 
NETWORKS 

by 

Lee,  Chan 

September  1983 

Thesis 

Advisor:                   R. 

Kevin  Wood 

Approved  for  public  release;  distribution  unlimited 


T214336 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whin  Dmtm  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  REPORT  NUMBER 


2.  GOVT  ACCESSION  NO 


3.  RECIPIENT'S  CATALOG  NUMBER 


4.     TITLE  (end  Subtitle) 


5.     TYPE   OF    REPORT    4    PERIOO   COVERED 


Efficient  Computation  of  the  K-terminal 
Reliability  of  Directed  Acyclic  Networks 


Master's  Thesis 
September  1983 


8.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AuTHORf*; 


Lee,  Chan 


S.  CONTRACT  OR  GRANT  NUMSEPrij 


t.  PERFORMING  ORGANIZATION  NAME  ANO  AOORE5S 

Naval  Postgraduate  School 
Monterey,  California  93943 


10.     PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  &    WORK  UNIT  NUMBERS 


II.     CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California  93943 


12.  REPORT  OATE 


13.  NUMBER  OF  PAGES 


47 


14.     MONITORING  AGENCY  NAME  *   AOORESS^/  dlfieront  Irom  Controlling  Office) 


15.     SECURITY   CLASS,   (of  thle  report) 


Unclassified 


15*.     DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


l«.     DISTRIBUTION   STATEMENT  (ol  thle  Xeport) 


Approved  for  public  release,  distribution  unlimited 


17.     DISTRIBUTION  STATEMENT  (ol  tho  eb, tract  entered  in  Block  20,  It  different  Irom  Roport) 


It.    SUPPLEMENTARY  NOTES 


It.     KEY  WORDS    Continue  on  tovtto  aide  II  nocooomty  mnd  Identify  by  block  number) 

Network  reliability,  Acyclic,  Directed  graph,  Vertex, 
Edge,  Reduction,  Decomposition,  Complexity 


20.     ABSTRACT  (Continue  on  rovoreo  etdo  tl  nocoeeewy  end  Identity  by  block  number) 


Several  topological  reduction  and  decomposition  techniques  are 
developed  to  decrease  the  complexity  of  computing  R„K(D),  che  source-to-K 

terminal  reliability  of  an  acyclic  directed  network  D  with  indeoendent 
component  failures.  R  K(D)  is  computed  in  O(IVIIEI)  time  wnen  D  can  be 

completely  reduced.  When  not  completely  reduced,  a  graph  D'  remains 


DO    ,  JAH  73    1473  EDITION  OF   1  NOV  «S  IS  OBSOLETE 

S/N  0102-  LF- 014-6601 


Unclassified 


1        SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dote  Enterec 


Uncassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGe  (Whmt  Dmtm  Ener»dj 


#20  -  ABSTRACT  -  (Continued) 

such  that  R  „(D)=M  R  mV   (DjR  ,,  (D0)  ...  R  .,  (Dm)  where  M  is  a  known  constant 

S  In  S  l\-i    I    S  l\o    ^-  S  In 

I       I  m 

and  the  D.  are   one  or  more  separable  components  of  D'.  A  simple  scheme, 

exponential  in  |V.-K.|  is  given  for  computing  R  ^(D).  When  1  V„.-K.|  become 

11  S  In  11 

too  large,  a  truncated  version  of  this  scheme  usually  gives  an  excellent 
lower  bound  on  R  ,  (D.)  and  thus  an  excellent  lower  bound  on  R,l/(D). 

S  In  •    1  S  In 

A  program  using  these  techniques  has  been  coded  in  FORTRAN  and  tested 
on  "complete"  acyclic  graphs  and  "street"  networks  with  up  to  100  vertices. 
Running  on  an  IBM  3033AP  under  FORTRAN  H  (Extended),  total  CPU  time  for 
computing  exact  reliability  is  less  than  3.5  seconds  when  IV-KI^.10  and 
IVI^IOO. 


S-N  0102-  LF-  014-  6601 

Unclassified 


j    SECURITY  CLASSIFICATION  OF  THIS  FAGefW>i#n  Dtm  Enfrmd) 


Approved  for  Public  Release;  Distributed  unlimited 

Efficient  Computation  of  the  K-Terminal  Reliability 
of  Directed  Acyclic  Networks 

by 


Lee.  Chan 

Major,  Republic  of  Korea  A.F. 

B.S..  Republic  of  Korea  A.F.  Academy,  1973 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
SeDtember  1983. 


37£ 
/ 


ABSTRACT 

Several  topological  reduction  and  decomposition  techniques  are  developed  to 
decrease  the  complexity  of  computing  R,x{D),  the  source-to-K-terminal  reliability 
of  an  acyclic  directed  network  D  with  independent  component  failures.  R5k(&)  is 
computed  in  0(  IVJI.E'I  )  time  when  D  can  be  completely  reduced.  When  not 
completely  reduced,  a  graph  D'  remains  such  that  RsK{D)  = 
M'RxxiD^RsiciDz)  ■  ■  ■  R2k  (An)  where  M  ls  a  known  constant  and  the  D-.  are  one 
or  more  separable  components  of  D' .  A  simple  scheme,  exponential  in  |  Mi  —  KL 
given  for  computing  Rsk(D)-  When  \Vt  —  A*  |  become  too  large,  a  truncated  version 
of  this  scheme  usually  gives  an  excellent  lower  bound  on  /?ajq(A)  an^  thus 
excellent  lower  bound  on  Rax(D). 

A  program  using  these  techniques  has  been  coded  in  FORTRAN  and  tested  on 
"complete"  acyclic  graphs  and  "street"  networks  with  up  to  ICO  vertices.  Running 
on  an  IBM  3033AP  under  FORTRAN  H  (Extended),  total  CPU  time  for  computing 
exact  reliability  is  less  than  3.5  seconds  when  i  V-Kl^lQ  and  j  7|<:0Q. 
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I.    INTRODUCTION 

Often,  in  the  design  and  production  of  complex  system  such  as  computer 
network,  good  system  reliability  is  a  critical  requirement.  Network  reliability  is 
rarely  treated  explicitly  and  quantitatively,  however,  largely  because  of  problems 
of  computational  complexity  and  because  few  efficient  computer  programs  exist 
for  computing  network  reliability.  Only  in  the  last  few  years  have  research  and 
development  in  this  area  provided  usable  results. 

We  wish,  as  efficiently  as  possible,  to  determine  the  network  reliability  for 
specific  class  of  networks.  The  purpose  of  this  thesis  is  to  develop  new  reduction 
and  decomposition  techniques  for  computing  network  reliability  and  to  show  how 
computational  complexity  can  be  decreased  by  using  these  techniques. 

The  network  model  which  is  used  in  this  thesis  may  be  thought  of  as  a 
communication  network  with  directed  edges,  allowing  only  one  way  communication 
from  tail  vertex  to  head  vertex.  Additionally,  we  impose  the  restriction  that  the 
network  have  no  directed  cycles.  Such  a  model  may  be  appropriate  in  hierarchical 
communication  networks  such  as  command  networks.  Also,  such  networks  may 
arise  as  subproblems  while  computing  the  reliability  of  general  directed  networks, 
Agrawal  [1982]. 

Networks  can  be  divided  into  various  classes,  e.g.,  directed,  undirected,  cyclic. 
acyclic,  etc.  A  large  amount  of  literature  exists  on  each  class.  (See  Agrawal 
[1982],  Ball  and  Provan  [1981],  Satyanarayana  [1982],  Hagstrom  [1980]  and 
Satyanarayana  and  Prabahkar  [1978]  for  directed  networks;  see  Ball  and  Provan 
[1981],  Johnson  [1982],  Satyanarayana  [1982],  Satyanarayana  and  Chang  [1981], 
Satyanaryana  and  Wood  [1982],  Valiant  [1978],   and  Wood  [1982]   for  undirected 
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networks;  see  Agrawal  [1982]  and  Satyanarayana  and  Prabhakar  [  1978]  for  cyclic 
versus  acyclic  directed  networks.)  The  specific  model  we  will  consider  is  the 
directed  acyclic  network,  denoted  D{V,E)  where  V  is  the  set  of  vertices  and  E  is 
the  set  of  edges.  Each  edge  is  represented  by  an  ordered  pair  of  vertices  e  =(u,v). 
The  ordering  implies  that  the  flow  of  communication  can  move  along  e  from  u  to  v 
only. 

Each  edge  e*  in  the  network  functions  with  probability  p^  and  fails  with 
probability  gt  =  l-Pi-  All  edge  failures  are  assumed  to  be  independent  events.  For 
simplicity  in  the  initial  discussion,  vertices  are  assumed  to  be  completely  reliable. 
However,  the  extension  to  unreliable  vertices  will  be  handled  briefly. 

Vertex  u  can  send  communication  to  vertex  v  in  D  if  there  is  a  directed  path 
of  working  edges  from  u  to  v.  One  of  the  vertices  s  is  designated  as  the  "source" 
or  the  "root".  With  respect  to  the  source,  there  are  three  different  measures  of 
reliability  which  are  usually  studied  in  directed  networks: 

(1)  The  source-to-terminal  (ST)  reliability,  Rst(D),  is  the  probability  that  s   can 

send  communication  to  a  specified  vertex  t  in  D. 

(2)  The  source-to-all-terminal  (SAT)  reliability,  Ray(D).  is  the  probability  that  «?  can 

send  communication  to  all  vertices  of  D. 

(3)  The  source-to-K-terminal  (SKT)  reliability,  R*g(D),  is  the  probability  that  s  can 

send  communication  to  every  vertex  in  a  specified  set  K  (these  vertices  are 

referred  to  as  "K-vertices")  with,  by  convention,  se/f. 
Of  course,  (l)  and  (2)  are  special  cases  of  (3).    We  consider  the  general  case  (3)  in 
this  thesis. 

Generally,  the  of  problem  computing  network  reliability  is  NP-hard.  Certain 
reduction  and  decomposition  schemes  exist  which  reduce  the  size  of  the  problem 
while  preserving  network  reliability,  e.g.,  Agrawal[l982],  Shogan[l973]  and  Wood 
[1982],   and,   for  special  classes  of  networks,   it  is  thereby  possible  to  compute 


reliability  in  polynomial  time.  In  general,  however,  the  solution  time  grows 
exponentially  in  the  size  of  the  network.  In  this  thesis,  we  develop  a  method  for 
computing  reliability  which,  in  general,  is  exponential  in  the  size  of  D ,  but  which  is 
polynomial  when  the  number  of  non-K-vertices.  i.e.,    |  /— K\ ,  is  fixed. 

The  rest  of  this  thesis  is  outlined  as  follows. 

Chapter  2  contains  a  necessary  definitions  and  graph  theoretic  concepts. 

Chapter  3  develops  reduction  and  decomposition  techniques  to  reduce  the  size 
of  a  network  in  order  to  decrease  computation  time.  Some  of  these  techniques  are 
new. 

Chapter  4  shows  how  to  compute  network  reliability  exactly  when  the  number 
of  non-K-vertices  is  moderate.  A  truncated  version  of  this  technique  is  used  to 
obtain  a  lower  bound  on  reliability  when  the  number  of  non-K-vertices  is  too  large. 
We  demonstrate  that  this  lower  bound  will  usually  be  "good." 

Chapter  5  describes  an  algorithm  to  compute  Rsx{D)  or  find  a  lower  bound  on 
Hsfr(D),  and  its  implementation.  The  program  is  coded  in  FORTRAN  with  one  main 
routine  and  18  subroutines  totaling  about  1600  lines  of  code. 

Chapter  6  gives  computational  experience  for  the  algorithm  run  on  an  IBM 
3033AP  system  /370  with  the  FORTRAN  H  (Extended)  compiler. 

Chapter  7  is  a  conclusion  and  gives  suggestions  for  further  research. 
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II.    DEFINITIONS  AND  NETWORK  RELIABILITY  THEORY 


In  this  chapter,  we  give  the  definitions  and  notation  used  throughout  this 
thesis,  and  introduce  basic  network  reliability  theory. 

A.      DEFINITIONS  AND  NOTATION 

A  directed  network  D—(V,E)  comprises  two  finite  sets:  V  is  the  set  of  vertices 
and  E  is  the  set  of  edges.  Each  edge  eeE  corresponds  to  an  ordered  pair  of 
vertices  in  the  directed  network.  That  is,  e={u,v)  where  u.veK.  The  vertex  u  is 
called  the  tail  vertex  of  edge  e  and  vertex  v  is  called  the  head  vertex  of  edge  s . 
For  any  two  vertices  u  and  v  in  D,  if  there  exists  an  edge  t=(u,v),  then  u  and  v 
are  said  to  be  adjacent  and  edge  e  is  incident  into  v  and  incident  out  of  of  u .  The 
indegree  of  a  vertex  v,  denoted  indeg(v),  is  the  number  of  edges  incident  into  v 
and  the  outdegree  of  a  vertex  v,  denoted  outdeg(v),  is  the  number  of  edges 
incident  out  of  v .  Two  vertices  are  connected  (or  "communicate")  if  there  exists  a 
sequence  of  vertices  and  edges  of  the  form  i>o.  (^o.^i).  vi.  {v \>vz) >■  ■■  >\vm-\>vm).  vm- 
This  sequence  is  a  path  of  length  m.   If  the  vQ=vm,  the  path  is  a  cycle. 

Consider  a  undirected  graph  G(V,E)  formed  by  ignoring  the  ordering  of  the 
edges  in  a  directed  graph  D{V,E).  G  is  called  the  underlying  graph  of  D.  D  is 
connected  if  its  underlying  graph  G  is  connected.  G  is  connected  if  there  is  an 
undirected  path  between  all  pairs  of  vertices.  The  vertex  connectivity  of  D  is  the 
minimum  number  of  vertices  which  must  be  deleted  from  G  (along  with  adjacent 
edges,  of  course)  in  order  to  disconnect  G  or  to  create  a  trivial  graph  with  only  one 
vertex.   D  is  separable  if  its  vertex  connectivity  is  one.   Otherwise  D  is 
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nonseparable.   In  a  separable  graph,  a  vertex  whose  deletion  disconnects  the  graph 
is  called  a  cut-verter. 

Merging  a  vertex  v  into  another  vertex  u>  means  removing  ail  edges  between  v 
and  w  and  forming  a  single  vertex  v'  by  coalescing  v  and  w  such  that  any  edge 
incident  into  (out  of)  v  or  w  is  now  incident  into  (out  of)  v'.  Merging  v  and  w  is 
sometimes  also  referred  to  as  "identifying"  v  with  w. 

The  SKT  reliability  of  a  directed  network  D,  denoted  RsK{D),  is  the  probability 
that  all  K-vertices  in  D  are  connected  by  path  of  working  edges  and  vertices  from 
source  vertex  s.  The  probability  that  edge  et  works  is  p+  =  l-Ji-  For  now,  we 
assume  that  all  vertices  work  with  probability  one. 

B.       NETWORK  THEORY  AND  SURVEY 
1.       XP-rnmplptpnpss 

We  briefly  discuss  the  complexity  of  the  SKT  reliability  problem  in 
terms  of  the  theory  of  NP-completeness  and  in  terms  of  practical  computation. 
More  detail  can  be  found  in  Garey  and  Johnson  [1979],  Ball  and  Provan[l98l]. 

Algorithms  to  solve  various  problems  can  be  broadly  classified  into  two 
categories,  polynomial- time  and  exponential-time  algorithms.  An  algorithm  is  a 
polynomial-time  algorithm  if  for  a  problem  of  size  n,  its  running  time  is  bounded 
by  a  polynomial  in  n.  Any  algorithm  that  is  not  a  polynomial-time  algorithm  is  an 
exponential-time  algorithm.  In  combinatorics,  a  decision  problem  P  is  said  to 
belong  to  the  class  NP  if  given  a  tentative  solution,  the  validity  of  the  solution  can 
be  checked  in  polynomial-time.  P  is  "NP-complete"  if  it  is  equivalent  to  the 
satisfiability  problem  of  a  conjunctive  normal  form  boolean  expression.  A  problem 
is  NP-hard  if  it  is  at  least  as  hard  as  an  NP-complete  problem. 
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A  problem  is  #P-complete  (number-P  complete)  if  its  solution  is  equivalent  to 
counting  (not  enumerating)  the  number  of  solutions  to  an  NP-comple':e  problem. 
Any  //P-complete  problem  is  NP-hard. 

2.  Pmhipm  Complexity 

Bail  and  Provan  [1981]  show  that  the  ST  reliability  problem  in  an  acyclic 
directed  network  is  a  #P-complete  problem.  This  implies  that  the  SKT  problem  is 
also  #P-complete  since  problem  instances  restricted  to  \K\=2  are  #P-complete. 
Even  finding  an  c-approximation  to  SKT  reliability  is  an  #P-complete  problem. 

In  practice,  these  complexity  results  imply  that  a  general  problem  of 
large  size  cannot  be  solved  analytically  in  a  reasonable  amount  of  time.  In  fact, 
the  time  taken  grows  exponentially  with  size  of  the  network.  As  will  be  seen, 
however,  practical  problems  of  fairly  large  size  can  be  handled  with  the  techniques 
developed  here.  Furthermore,  even  larger  problems  can  be  handled  if  one  is 
willing  to  accept  a  "good"  lower  bound  on  exact  reliability. 

3.  Seduction 

Many  reduction  schemes  have  been  developed  to  help  solve  various 
network  reliability  problems.  These  reductions  reduce  the  size  of  the  network  in 
polynomial  time  while  preserving  reliability,  i.e.,  they  create  a  new,  smaller, 
network  G'  from  the  original  network  G  [G  may  directed  or  undirected  and  any 
number  of  definitions  of  reliability  may  be  used.)  such  that  R{G)~M  K(G')  where  M 
is  a  known  constant  factor.  In  this  way,  a  polynomial  amount  cf  work  is  expended 
in  order  to  save,  it  is  hoped,  an  exponential  amount  of  work  during  later,  more 
general  computations.  In  some  special,  but  often  practical  cases,  it  is  possible  to 
completely  reduce  a  network  and  effectively  compute  its  reliability  in  polynomial 
time. 

The  reductions  describe  below  give  a  flavor  for  the  different  types  of 
reductions  which  have  been  developed.  Parallel  and  series  reductions  are  well- 
known  schemes  so  they  will  not  be  discussed.    The  reductions  described  assume 
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that  vertices  are  completely  reliable.  However,  most  can  be  extended  to  handle 
unreliable  vertices. 

(1)  Polygon-To-Chain  Reductions:  A  chain  C  in  an  undirected  graph  G  is  defined  as 
a  alternating  sequence  of  vertices  and  edges  beginning  and  ending  with  vertices 
such  that  ail  internal  vertices  have  degree  2.  Two  parallel  chains  form  a  polygon 
and  this  polygon  can  be  replaced  by  a  chain  in  a  reduction  that  preserves 
reliability  in  the  K-terminal  reliability  problem,  the  analog  to  the  SKT  reliability 
problem  in  an  undirected  network.  This  reduction  and  its  implications  are 
developed  by  Yfood  [1982]. 

(2)  2-Neighbor  Vertex  Reduction:  For  the  SKT  reliability  problem,  any  vertex  u 
that  has  exactly  two  neighbors  v  and  w  can  be  eliminated  from  the  network  if 
ur/  K  or  \u,v ,w\zK  and  u*s  . 

(3)  Parallel-Chain  Reduction:  Suppose  Eu  is  the  set  of  all  edges  in  a  directed 
network  D  whose  endpoints  are  ty  and  Vj.  In  D,  a  chain  C  is  an  alternating 
sequence  of  vertices  and  non-empty  sets  of  edges  where  the  internal  vertices  are 
2-neighbor  vertices  in  D.  Two  chains  C,  and  C2  with  identical  endpoints  are 
parallel  chains.  Parallel  chains  can  be  reduced  to  a  single  chain  in  the  SKT 
reliability  problem.   Reductions  (2)  and  (3)  are  developed  by  Agrawal  [1932]. 

4.       Decomposition 

Sometimes  it  is  possible  to  divide  a  given  graph  into  a  number  of 
subgraphs,  such  that  each  subgraph  can  be  analyzed  separately  and  the  results 
combined  to  obtain  the  reliability  of  the  whole  graph.  We  briefly  describe  some  of 
these  methods  based  on  vertex  connectivity. 

(1)  2-Connected  Graph  Decomposition  Let  G  be  a  2-connected  graph  with  a  pair  of 
vertices  \u,v\  whose  deletion  disconnects  the  graph.  These  vertices  are  called  a 
"separating  pair."  Tins  separating  pair  partitions  G  into  two  subgraphs  Gl  and  G2 
such  that  Gj^Gj  =G,  Gi^Gz  -\u,v\,  E^^  and  Ez*$-  It  has  been  shown  by  Wcod 
[1982]  that,  in  the  K-terminal  reliability  problem,  each  subgraph  defines  a  chain 

12 


through  a  set  of  reliability  computations  on  the  subgraphs.  The  subgraphs  may  be 
replaced  by  the  chains  so  as  to  preserve  reliability.  Agrawal  [1982]  has  proven  the 
analogous  decomposition  scheme  for  the  SKT  reliability  problem  on  directed 
networks.  (2)  Strong-connected  Decomposition:  A  strong  component  of  a  directed 
graph  D  is  a  maximal  subgraph  of  D  such  that  there  exists  a  directed  path  from 
every  vertex  to  every  other  vertex  in  the  subgraph.  w"hen  computing  Rsx(D),  each 
strong  component  can  be  considered  separately.  Computing  each  component 
#»tf(A)  and  multipling  /?«r(A)  is  whole  graph  reliability  Rsk{&)-  This 
decomposition  has  been  developed  by  Shogan  [1973]  and  Agrawal  [1982]. 

Other  decomposition  schemes  exist.  Decomposition  in  3-connected 
graphs  is  often  possible,  Rosenthal  [1974].  Decomposition  in  1-connected  graphs 
will  be  discussed  in  chapter  3. 
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III.    REDUCTION  AND  DECOMPOSITION 

Any  algorithm  to  solve  a  general  network  reliability  problem  has  a  running 
time  which  is  exponential  in  the  size  of  problem.  If  by  some  method,  the  size  of  the 
problem  can  be  reduced,  immense  computational  savings  can  be  made.  Therefore, 
in  this  chapter,  we  introduce  reduction  and  decomposition  schemes  which  will 
often  reduce  the  size  of  the  problem  while  always  preserving  reliability.  Some  of 
these  techniques  are  specific  to  the  SKT  reliability  problem  in  acyclic  directed 
networks  while  others  are  of  more  general  applicability. 

Reduction  schemes  reduce  a  parameter  of  the  size  of  the  problem,  usually 
edges  or  vertices,  in  order  to  decrease  computational  complexity  while  preserving 
network  reliability.  On  the  other  hand,  decomposition  schemes  split  the  whole 
graph  into  subgraphs,  computes  subgraphs  reliabilities  separately  and  then 
combines  these  reliabilities  in  some  way  to  compute  overall  network  reliability. 

A.      REDUCTIONS 

In  this  section,  we  discuss  various  reduction  schemes  which  reduce  the  size  of 
D  and,  consequently,  reduce  the  complexity  of  computing  Rsk{D).  Under  any 
reduction,  reliability  of  the  network  remains  invariant  up  to  a  known,  constant 
factor,  i.e.,  BsK(D)  =  M*RaiC{D')  where  M  is  the  known  factor,  FT  is  the  new  set  of 
K-vertices,  and  D'  is  the  reduced  network.  In  the  following  illustrations,  all  vertices 
are  represented  as  circles  with  K-vertices  being  darkened  or  shaded. 

The  parallel  and  series  reductions  are  standard,  i.e.,  well-known  and  widely 
used.   The  three  reduction  schemes, 
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"indegree-one  K-vertex  contraction".  "non-K-vertex  deletion"  and  "neck  vertex  to 
K-vertex  reduction"  are  newly  developed  reductions. 

1.  Par^llpi  FHgpq  RpHnrMnn 

In  a  directed  network,  two  parallel  edges  ea=(u,v)  and  eb=(u,v)  may 
be  replaced  a  single  edge  ec=(u,v)  while  preserving  reliability  if  pc  -  pa+Pb-Pa?b 
andJ/=l. 

2.  SpHp<?  Fl^cpc;  PpHnrtinn 

In  a  directed  network,  a  pair  of  edges  ea={u,v)  and  eb={v,w)  such  that 
indeg(z;)=outdeg(T;)  =  l  are  series  edges.  If  v&  K,  edges  e0  and  e6  can  be  replaced 
by  single  edge  ec  =  {u,w)  while  preserving  reliability  if  pc  =  papb  and  if  si. 

3.  TnHpgrpp-OnP  K-wrfpy  fnnt-rarnnn 

Let  ei  =  (u,v)  be  an  edge  in  D  such  that  ve.K  and  indeg(v)  =  l.  In  the  5KT 
reliability  prcblem.  communication  from  s  to  vertex  vzK  can  only  only  occur  via 
edge  ei  so  this  edge  must  work  if  the  network  is  to  work.  Consequently,  edge  e* 
can  be  contracted  such  that  R+(D)  =  M-R3k(D')  and  M-Pi 

where  D'  -  \  V  —  u  —  v  r  w  ,E  -  et } 
w  =  u\jv 


K  - 


K—v  —a  +w ,  if  u  €.K. 


The  validity  of  this  reduction  is  easily  shown.  Let  et(u,i/)  be  the  edge 
described  above  in  the  acyclic  network  D.  Let  F\  denote  the  event  that  e^  is 
working,  and  let  Fi  denote  the  complementary  event.  Since  R3x(D)  is  just  a 
probability,  the  rules  of  conditional  probability  can  be  applied  to  obtaun 

•      Rsjc(D)=piRsK(D\Fi)  +  qiRsjC(D\Fi) 

=  piR,K{D\Fi)    since  R^iDlF^O 
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since  no  new  paths  are  created  by  contraction. 

4.  Mnn-fC-VprfpY  Dplprinn 

In  the  SKT  reliability  problem,  an  outdegree  zero  non-K-vertex  can  be 
deleted  from  the  network  without  changing  network  reliability  since  such  a  vertex 
is  irrelevant.   In  this  case  Rtx(D)  -  M  Rsk(D-v)  and  Afsl. 

5.  Nprk  Vprtpy  tr>  K-vprtPY  Rprhirtinn 

Suppose  u  and  v  are  two  vertices  in  D  such  that  uzK,  v3  K  and  all 
paths  from  s  to  u  include  v.  v  is  called  a  neck  vertex.  It  may  or  may  not  be  a 
cut-vertex.  Since  s  must  communicate  with  v  before  it  communicates  with  u, 
RsK(D)  =  M  /?s(^>y)(.D)  where  M-\.  Thus,  the  neck  vertex  to  K-vertex  reduction 
simply  changes  K  to  K+v .  See  Figures  3.1a  and  3.ib.  (All  figures  are  given  in 
Appendix  A.) 

The  neck  vertex  to  K-vertex  is  termed  a  "reduction"  because  it  reduces 
the  parameter  |  V-K\  in  D.  As  vail  be  seen  in  chapter  4.  the  general  algorithm  for 
computing  Rsk(D)  runs  in  time  proportional  to  2^V~K^  so  this  is  a  reduction 
according  to  our  definition.  Furthermore,  the  neck  vertex  to  K-vertex  reduction 
will  often  allow  other  reductions  or  decomposition  to  take  place. 

6.  Qthpr  ffpHnrtinnq 

Numerous  other  reductions  exist  as  exemplified  in  chapter  2,  section 
B3.  However,  none  of  these  reductions  has  been  implemented  because  of  their 
complexity,  and  therefore,  no  other  reductions  will  be  described  here. 

B.     DECOMPOSITION 

Sometimes  it  is  possible  to  decompose  a  given  network  D  into  a  number  of 
subnetworks.  Each  subnetwork  can  then  be  analyzed  separately  and  the  results 
combined  in  some  way  to  obtain  overall  network  reliability. 
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Decomposition  can  result  in  substantial  savings  because,  in  general,  the  time 
to  analyze  the  whole  system  is  exponential  in  size  of  the  system.  Effort  needed  to 
compute  R,k(D)  using  decomposition  is  the  sum  of  efforts  needed  for  each 
subproblem  plus  the  effort,  normally  polynomial  in  the  size  of  D ,  in  performing  the 
decomposition  and  recombining  the  result.  Without  decomposition  the  necessary 
work  is  roughly  the  product  of  the  effort  needed  to  solve  the  individual  problems. 

1.  Hnt-Vprt-gv  Dprnrnpngitinn 

A  connected,  undirected,  network  G{  V,E)  has  a  cut-vertex  v  if  there 
exist  components  (subgraphs)  Cx  and  C2  such  that  ClDC2=v  and  C1^C2=  G  .  Any 
communication  from  one  component  to  the  other  must  pass  through  v.  In  this 
case  we  aiso  say  that  v  separates  Cx  and  C2-  A  network  which  has  a  cut-vertex  is 
"separable"  and  one  which  has  none  is  called  "nonseparabie."  A  directed  network  D 
is  separable  if  its  underlying  graph  G  is  separable. 

Let  D  have  nonseparabie  components  Cl,C2,  ■  ■  .  ,  C»-  Any  component  Q 
which  has  only  one  cut-vertex  is  called  a  pendant  component.  If  a  pendant 
component,  Q,  contains  no  K-vertices  except  possibly  its  cut -vertex,  then  it  is 
completely  irrelevant  with  respect  to  reliability  since  none  of  its  edges  are  in  any 
path  from  the  source  to  a  K-vertex.  Thus,  this  component,  excluding  its  cut-vertex, 
can  be  deleted  from  the  network  while  preserving  network  reliability  in  the  SKT 
reliability  problem.  This  type  of  deletion  may  be  recursively  applied  until  only 
relevant  components  remain. 

Cut -vertices  are  not  efficiently  found  in  a  directed  network  if  each  edge 
is  represented  in  one  direction,  e.g.,  if  all  edges  incident  into  each  vertex  are 
stored  contiguously.  However,  if  edges  are  stored  both  by  tail  vertex  and  head 
vertex,  then  this  is  exactly  the  representation  of  the  underlying  undirected 
network.  Thus,  an  efficient  algorithm  for  finding  cut -vertices  in  an  undirected 
graph  may  be  used  to  find  cut-vertices  in  D.  A  slight  modification  to  the 
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DFS  (depth-first  search)  algorithm  due  to  Tarjan  [1972]  can  find  ail  cut-vertices  in 
Q(\&\)  time. 

Once  all  cut-vertices  have  been  found  and  all  irrelevant  components 
removed,  any  cut -vertices  remaining  which  are  not  K-vertices  are  neck  vertices 
and  may  be  added  to  the  set  K.  It  is  fairly  obvious  that  each  component  must  work 
by  itself  if  D  is  to  work.  Since  edge  reliabilities  are  independent  it  follows  that 
RSK{D)-Rix{Dl)Rzx{D2)  •  •■  Rgxi^m)  where  the  Di  are  the  nonseparabie, 
relevant  components  of  D.  This  is  identically  equal  to  a  network  which  has  the 
same  components  as  D  but  whose  only  cut-vertex  is  the  source.  A  more  rigorous 
demonstration  of  this  identity  can  be  made  by  applying  the  decomposition  scheme 
described  next,  "moving  edges  to  source,"  to  all  edges  incident  out  of  cut-vertices. 
Figure  3.2  demonstrates  the  results  of  cut-vertex  decomposition  applied  to  a 
network  D  with  five  components,  two  of  which  are  irrelevant.  In  this  case, 
Rsx(D)  =  RsK(Dl).RsK  (D2)-RsK(D5). 

2.  Moving  GdftCS  LQ  ^,n\-\rr>o 

In  this  section,  we  discuss  a  technique  which  does  not  fit  neatly  under 
the  rubric  of  decomposition  or  reduction  since  it  only  moves  edges  from  one  place 
to  another  in  the  network.  However,  by  moving  certain  edge  tails  to  the  source  s, 
new  parallel  reductions,  series  reductions  or  contractions  may  be  made  possible 
along  with  cut-vertex  decompositions.  Consider  an  edge  e-{u,v)  such  that  uzK. 
Edge  e  can  be  made  incident  out  of  any  other  K-vertex  il'zK  without  affecting 
reliability  as  long  as  no  cycles  are  created.  A  simple  conditioning  argument  can  be 
-used  to  show  th±s,  Agrawal  [1982].  The  source  is  the  obvious  candidate  to  receive 
the  edge  tails  since  it  has  no  edges  incident  into  it  and  therefore  no  cycles  can 
possibly  be  created.  Figure  3.3  shows  an  example  where  this  operation  creates  a 
cut-vertex  and  allows  a  parallel  reduction. 
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IV.    COMPUTING  ag  BUittUDL  B  AN  ACYCLIC  DIRECTED  NETTORK 

This  chapter  gives  a  general  method  for  computing  SKT  reliabiLty  in  an  acyclic 
directed  network.  The  method  is  based  on  repeated  computation  of  SAT  reliability 
of  certain  subgraphs  of  D  and,  consequently,  we  first  discuss  SAT  reliability.  After 
deriving  a  method  for  computing  SKT  reliability,  a  simple  lower  bound  becomes 
apparent  and  several  heuristic  arguments  are  given  indicating  that  the  bound's 
accuracy  should  be  good.  Actual  computational  examples  are  reserved  for  the 
next  chapter. 


A.      SOURCE-TO-ALL  TERMLNAL  RELIABILITY 

Let  D{V,E)  be  an  acyclic  directed  network  and  let  Jt=l-JK  be  the  failure 
probability  for  any  edge  2iEE  and  let  Ej  be  the  set  of  edges  incident  into  Vj.  Also, 
let  the  vertices  be  numbered  such  that  [1/1,1/2,  •  ■  •  ,vn  j  is  an  an  acyclic  ordering  of 
V.   Then.  SAT  reliability  may  be  expressed  as 

Rsy{D)  =  P{  s  communicates  with  all  v  e  V^vn  (4.1) 

and  at  least  one  edge  into  vn  works  ) 
since  there  no  edges  incident  out  of  vn  by  the  acyclic  ordering 

by  independence 


=  Rs{V-vj{D-njn) 


1-  n  n 


I"  H  7; 


nf:-n  g, 
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Another  way  to  see  this  is  to  use  the  fundamental  topological  property  that 
RsY{D)  is  unchanged  if  any  edge  (u,v)  with  ueiT  is  replaced  by  an  edge  {w,v)  with 
the  same  reliability  and  weK  has  a  lower  number  than  u  in  an  acyclic  ordering  of 
D  Consequently,  in  the  SAT  problem,  we  may  move  ail  edges  so  that  they  come 
out  of  the  source  just  like  the  "moving  edge  to  source"  operation  in  section  3.2. 
Equation  (4-.1)  is  then  trivially  true  for  a  network  where  all  edges  come  out  of  the 
source.  (See  Figure  4.1).  The  complexity  of  computing  Rsv{D)  Is  easily  seen  to  be 
0(  \E  |)  if  all  edges  are  with  ccmmon  head  vertices  are  stored  contiguously. 

B.      SOURCE-TO-K-TERMINAL  RELIABILITY 

Now  we  consider  the  SKT  reliability  problem  in  D.  Let  D{V,E)  be  a  acyclic 
directed  network  with  KqV.   By  simple  state-space  partitioning, 

ftsk  (-C)  =  P{  s  communicates  with  all  v  s  V) 

+    S  P(s  communicates  with  all  vz  V-v^  but  not  with  vj 

"**      S     P(  s  communicates  with  all  v^V-v^—Vj  but  not  with  Vi  or  Vj) 
+  .  .  . 

+  P{s  communicates  with  all  v  eK  but  not  with  any  v  e  V-K) 

But, 

P{  s  communicates  with  all  vzV)  =  RsV{D) 

P{  s  communicates  with  all  vsV—Vi  but  not  withi^   )  =  Rs{v-v j(^)-  P(  ^°  edge  into 

Vi  works  ) 
P(  s  communicates  with  all  ye7-ui-i;J)=i?,(^t,1  _..)(/?).  P{  No  edge  into  U|  or  Vj 

works  except  an  edge  e=(t»i,i;J)  or  e  =(1^,1^)) 
etc., 
In  general,  then. 
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RsK{D)  =  RsV(D)  +    EJ[.n  <h]*.AD  -^)|  (4.2) 


+       .     . 


where  E*  is  the  set  of  edges  going  into  vit  Eq  is  the  set  of  all  edges  going  into  v<  or 
Vj  except  edges  going  between  vt  and  Vj,  Eqt  is  the  set  of  all  edges  going  into  vit  Vj 
or  Vj.  except  edges  going  betv/een  pairs  of  vertices  in  \vitVj.vk],  etc. 

The  complexity  of  computing  R^iD)  is,  via  the  above  formula,  exponential  :n 
the  number  of  non-K-vertices.   The  total  number  of  additive  terms  in  equation  (4.2) 

is  2a  0!V~K^  =  &V~K][  where  C?  is  the  number  of  ways  to  choose  i  elements  out  of  n 

i=0 

elements.  Computing  SKT  reliability  thus  requires  Q(\E\2^V~K])  time  using 
equation  (4.2).  Of  course,  if  the  number  of  non-K-vertices  is  fixed  as  \E\  changes, 
computation  is  effectively  0(  |  E  \ ). 

C.     A  LOVER  BOUND  ON  SKT  RELIABILITY 

Since  computation  increases  exponentially  in  the  worst  case,  we  may  be 
satisfied  with  an  approximation  to  SKT  reliability  which  requires  less  time  to 
compute  than  exact  reliability.  Since  all  terms  in  equation  (4.2)  are  positive,  any 
subset  of  these  terms  provides  an  approximation  to  SKT  reliability  which  is  a  lower 
bound.  Also,  the  product  formed  by  lower  bounds  on  independent  components  of  a 
network  yields  a  lower  bound  on  overall  network  reliability.  The  subset  of  terms  we 
choose  to  form  ohis  lower,  bound  consists  of  terms  enumerated  in  the  order  given  in 
equation  (4.2)  until  some  specified  computational  limit  is  reached.  The  motivation 
for  this  scheme  being  fairly  accurate  is  given  below. 
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Suppose  Pi  =  p  for  all  eteE  and  p  is  close  to  1.    Consider  the  first  few  terms  of 
Equation  (4.2). 


RsV(D)=   nfl-g1*'1]    -    1.0 


IAI+I&I-I 


-1 


VjeV— 1/{— i/y 


I-,1*' 


ItfJH^I+l 


r  nji|%^,(H-*i)^,a,nil 

Thus,  the  contribution  of  the  terms  to  the  total  reliability  sum  diminishes 
rapidly  unless  there  are  a  significant  number  of  edges  connecting  non-K-vertices. 

Another  reason  that  the  lower  bound  should  often  be  good  is  that  terms  with 
many  vertices  deleted  may  be  identically  zero.  For  example,  consider  an  ST 
reliability  problem  where  the  shortest  path  from  s  to  t  has  L  edges  in  it.  This  path 
must  contain  l-l  non-K-vertices  and  so  Ra^y^Ui^v      ^Vk){D-/vi-Vj-...-^vk)  =  0  when 

|  V-\Vi,Vj vts,t\  |<i— 1.    Of  course,  other  terms  may  be  identically  zero,  too. 

For  example,  in  any  K-terminal  problem,  the  term  corresponding  to  D-^v^-njj  will 
be  identically  zero  if  all  paths  from  s  to  some  K-vertex  v  must  include  either  Vi  or 

vr 

Computational  experience  in  Chapter  V  lends  credence  to  the  accuracy  of  this 
lower  bound. 
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V    ALGORITHM  AND  DATA  STRUCTURE 


The  algorithm  presented  here  for  reliability  evaluation  of  directed  acyclic 
networks  is  based  on  the  theory  developed  in  chapter  3  and  chapter  4.  The 
program,  is  intended  for  general  use  so  it  was  coded  in  the  widely  available 
language,  FORTRAN.  This  language  is  available  on  most  computer  systems  and  is 
still  one  of  the  most  popular  languages  among  operations  research  analysts. 

A.      ALGORITHM 

The  objective  in  this  section  is  to  develop  an  efficient  algorithm  for  computing 
SKT  reliability  in  a  directed  acyclic  network.  Initially,  polynomial-time  reductions 
and  decomposition  are  used  to  reduce  network  parameters  and  decompose  the  the 
network  into  subnetworks,  if  possible.  If  reliability  is  not  completely  computed  via 
reliability-preserving  reductions,  the  general  SKT  reliability  expression  is 
computed  as  described  in  chapter  4.  If  a  specified  CPU  time  limit  is  reached  for 
any  component,  the  computation  is  truncated  and  a  lower  bound  on  reliability  is 
produced. 

There  is  a  one  main  algorithm  and  numerous  subroutine  algorithms  in  the 
program.  We  discuss  the  main  program  and  several  important  subroutina 
algorithms  for  understanding.  The  structure  of  the  main  program  is  based  on  a 
step  by  step  approach  which  avoids  redundant  work.  Comments  will  be  given  in 
curly  brackets. 
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Main 

Input:  A  directed  acyclic  network  D  with  vertex  set  V,  \  V\>2,  edge  set  E,  -E  ^2, 
and  set  KqV,  \K\>2.   Edge  reliability  pt  for  each  edge  et  in  E. 

Output:  Rsk(P)>  exact  reliability,  if  CPU  time  less  than  specified  amount,  otherwise 
a  lower  bound. 

(1)  ^Initialize}  M  =  1.0,  bound=. false. 

(2)  Read  all  edge  data  and  create  network  data  structure 

(3)  \ Change  neck  vertices  to  K-vertices]  K  =   Neck  (D,K) 

(4)  ^Perform  reductions  for  whole  network]   Reduct  (D,K,M) 

(5)  If  remaining  vertices  |  V\  £  1,  go  to  end. 

(6)  \DltD2 A]  =  BecompCO) 

(7)  For  each  subnetwork  A 

(a)  Unitialize]  lowbnd=. false.,  set  time  =  0. 

(b)  £Move  edge  tails  to  source]  For  each  edge  e=(vitVj)  with  v^K 
$move  Vi  to  .s  ]   Move(e) 

(c)  Construct  the  data  structure  for  subnetwork  A 

(d)  ^Reduce  subnetwork]   Redact  (DitKi,M) 

(e)  [  Initialize  ]  sumpro=0 

(f)  \  Find  SAT  reliability  ]  sumpro  =   Rsa(A). 

(g)  Count  the  remaining  non-K-vertices  for  Dit  i.e.,  if  vertex  number 
>  0,  count =count-rl 

(h)  If  counts  l,  \  find  the  terms  in  J]  (ngi)/?sa(Z?^i -vi)\  ] 

1)  [Find  the  next  combination  of  non-K-vertices] 

\vWvh wiJ  =   CombGV0 

2)  \  FindHgt   for  Eili%...k  ]   St  =  Multi^.^ vk) 

3)  sumpro=sumpro  +  M'-  Rsa(Z?— u;i-v;-  — -VjJ 

4)  If  elapsed  CPU  time  >  time  limit  then  lowbnd=.true.  and 
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bound =.  true. 

5)    If    lowbnd=. false,    and    all    combination    have    not    been 
enumerated  then  go  to  (8.h.  i) 
(i)  M-M  x  sumpro 
(8)  If  bound=. false,  print  "Exact  reliability  is"  M,  otherwise  print  "Lower 
bound  on  reliability  is"  M 
End  Main 

function  Neck  (D,K) 

Input:  D,K 

Output:  New  set  of  K-vertices 

(  This  routine  finds  all  neck  vertices  and  changes  them  to  K-vertices  ) 

(1)  Put  all  v^K  into  Qn 

(2)*  =  |A1 

(3)  While  <&*$ 

(a)  Remove  v  from  Qn 

(b)  Search  all  K-vertices  in  D-Vi  using  Depth-First  Search 

(c)  k'=  number  of  K-vertices  reached 
(c)  llk'Kk  then  K=K\J \v\ 

(3)  Return  (K). 
End  of  Neck 

subroutine  lted\ict(D,K,M) 

Input:  D,K,M 

Output:  Reduced  network  D  with  modified  K  and  M 

(    This    routine     performs    all    non-K-vertex    reductions,     series -parallel 

reductions  and  indegree-one  K-vertex  contractions  ) 

(1)  Construct  Qd  -  \  veV-K  \  outdeg(?;)=0  \ 
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(2)  While  Qd*$ 

(a)  For  each  vertex  Vj  incident  to  ut,  if  outdeg(uj)  =  l  and  i/ycK— K 
then  add  Vj  to  Qd . 

(b)  Nonkre  ( i/t) 

(3)  Construct  Qa  =  \  v^V-K  |  indeg(i;i)=  outdeg(i»t)=l  \  and  Qc  =  \  v^K  \ 
indeg(iO  =  l  1 

(4)  For  each^t 

(a)  Check  all  outgoing   edges,   if  e  {vi.Vj) ,e(vi,Vj)   are   found  then 
Parall(e) 

(b)  If  parallel  reductions  creates  new  series  or  contraction  vertices, 
put  these  into  Qs    or   Qc,  respectively. 

(5)  While  Qs  *$  and  Qc  *0 

(a)  If  &*0   Series(  Q,(i)) 

1)  If  series  reductions  create  new  parallel  edges,  do  parallel 
reduction.    Parall(e) 

2)  If   parallel  reductions   create   new   series   or   contraction 
vertices,  put  these  into  Qs  or  Qc ,  respectively. 

(b)  If  Qc  *0   Contra(  <i(i)) 

1)  If  the  contraction  creates  a  new  contraction  vertex,    put 
these  into  Qc 

2)  If  the  contraction  creates  new  parallel  edges,    do  parallel 
reduction  Parail(2 ) 

3)  If   parallel  reductions   create   new   series   or   contraction 
vertices,  put  these  into  Q3  or  Qs,  respectively 

(6)  Return  {D.K.M) 


Etid  Reduct 
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subroutine  Decomp(Z?) 

Input:  D 

Output:  Subnetworks  D^.Dz Dm 

(1)  While  not  ail  cut-vertices  have  been  found 

(a)  Using  DFS,  find  pendant  component  D^  with  cut-vertex  v< 

(b)  If  Di  has  no  K-vertices  except  possibly  v<  then  delete  Di 

(2)  Return  (Dl%  D2 Dm) 

End  Decorap 

After  reading  in  the  data  and  creating  the  neiwork  data  structure,  the 
algorithm  begins  by  changing  any  neck  vertices  to  K-vertices  by  a  call  to  Neck  at 
step  (3).  Next,  a  cail  to  Reduct  is  made  at  step  (4).  Reduct  first  removes  all 
outdegree  zero  non-K-vertices.  Through  the  queue  mechanism,  it  is  ensured  that 
all  initial  outdegree  zero  non-K-vertices  are  deleted  along  with  any  which  are 
created  as  the  reduction  proceeds.  A  similar  method  is  used  for  the  other 
reductions  schemes. 

The  other  reductions,  indegree-one-K-vertex  contractions  and  parallel  and 
series  reductions,  may  recursively  create  new  reducible  vertices  or  edges.  The 
newly  created  reducible  vertices  are  put  into  a  queue  instead  of  immediately  being 
reduced  to  avoid  redundant  work.  Parallel  edges  are  reduced  immediately, 
however.  The  network  of  Figure  5.1  is  completely  reduced  by  reductions  alone  and 
illustrates  how  a  reduction  can  create  other  reducible  edges  or  vertices.  A  call  to 
Reduct  would  result  in  the  set  of  reductions  given  below.  This  network  would  be 
completely  reduced  and  its  reliability  completely  calculated  by  these  reductions. 
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(1)  Reduce  parallel  edges 

(2)  contract  vertex 

(3)  series  vertex 

(4)  contract  vertex 

(5)  reduce  parallel  edges 

(6)  series  vertex 

(6)  series  vertex 

(7)  reduce  parallel  edges 

(8)  series  vertex 

(9)  reduce  parallel  edges 

(10)  contract  vertex 

(11)  reduce  parallel  edges 

(12)  contract  vertex 

(13)  reduce  parallel  edges 

(14)  contract  vertex 

After  reductions  are  performed,  a  network  will  decompose  into 
subcomponents  if  it  contains  any  cut -vertices.  This  is  done  at  step  (6)  of  Main  by  a 
call  to  Decomp.  The  rest  of  the  calculations  are  performed  iteratively  for  each 
component  D\  under  step  (7).  First,  all  edges  with  tail-vertices  in  K  are  moved  to 
the  source.  This  may  produce  additional  reductions,  so  another  call  to  Reduct  is 
made  for  the  subcomponent.  Then  the  terms  corresponding  to  Equation  4.2  (for 
that  component)  are  produced  until  all  such  terms  have  been  enumerated  or  the 
CPU  time  limit  is  reached.  Note  that  this  CPU  limit  is  based  on  the  CPU  time  used 
in  calculating  subcomponent  reliability  and  is  not  based  on  total  CPU  time. 

If  the  CPU  time  limit  is  reached  for  any  of  the  subcomponents,  the  final  value 
obtained  by  the  algorithm  is  a  lower  bound  on  reliability.  Otherwise,  the  value 
obtained  is  exact  reliability.    Choosing  a  CPU  time  limit  will  depend  on  computer 
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facilities  available  and  the  network  under  analysis.  Figure  6.5  shows  the  reliability 
increment  as  a  function  of  CPU  time.  The  lower  bound  obtained  using  a  reasonable 
CPU  time  limit  will  not  differ  much  from  the  exact  reliability  for  this  large  scale 
network  example 

B.     DATA  STRUCTURE 

There  a  many  ways  to  represent  a  network  on  a  digital  computer.  The  most 
straightforward  is  to  use  an  adjacency  matrix:  If  vertices  ty  and  Vj  are  adjacent, 
M(i,j)=l.  and  M(i,j)=0  otherwise.  Matrices  of  the  above  type  are,  in  practice, 
extremely  sparse  and  this  representation  is  very  inefficient  in  terms  of  space.  It  is 
often  inefficient  in  terms  of  execution  time,  too,  particularly  when  one  is 
interested  in  operations  that  require  retrieving  information  from  all  vertices 
adjacent  to  a  particular  vertex,  e.g.,  Neck,  Decomp,  Reduct,  and  Rsa.  In  this 
situation,  since  most  of  the  elements  of  M(i,j)  are  null,  the  algorithm  will  spend  a 
great  amount  of  time  retrieving  and  comparing  zero  values.  For  a  static  network, 
an  efficient  representation  is  a  packed  matrix.  In  our  algorithm,  the  network  is 
dynamic,  since  it  is  reduced  and  decomposed  many  times  during  execution. 
Therefore,  a  multi-linked  data  structure  (multiple  linked  lists  with  links  to  other 
data  structures  such  as  arrays)  is  used  to  represent  the  network.  Every  vertex  has 
a  linked  list  of  adjacent  vertices,  which,  besides  informing  which  vertices  are 
adjacent  to  it,  also  tells  whether  or  not  the  vertex  belongs  to  set  K.  Two  such  lists 
are  kept  for  each  vertex,  an  "adjacent  out"  list  and  an  "adjacent  into"  list.  Of 
course,  each  adjacency  entry  corresponds  to  an  edge.  So,  in  addition,  there  are 
pointers  indicating  the  addresses  where  information  about  edges  is  kept. 

Figure  5.3  illustrates  this  data  structure.  The  vertex  vx  only  can  send 
communication  to  vertices  v2,  and  1/3,  and  vertex  v2  can  receive  the 
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communication  from  vertex  vl  and  send  the  communication  to  vertices  v3  and  vA, 
and  so  on   This  structure  is  a  very  efficient  data  structure  for  dynamic  networks. 

C.     ALGORITHMIC  COMPLEXITY 

It  is  important  to  determine  the  complexity  of  the  various  parts  cf  the 
algorithm  so  that  overall  complexity  can  be  understood.  In  this  section,  we  analyze 
the  worst-case  complexity  of  the  subroutines,  as  written,  and  thus.  And  the 
complexity  of  the  overall  algorithm.  Not  all  of  the  subroutines  have  been  written 
as  efficiently  as  possible  because  cf  the  difficulty  in  programming  such  routines 
and  because  of  limited  time.  Furthermore,  for  problems  of  moderate  size,  it  is 
unlikely  that  much  efficiency  would  be  gained  using  more  sophisticated  routines, 
since,  in  practice,  most  of  the  routines  written  seem  to  operate  mere  efficiently 
than  the  worst-case  analysis  indicates. 

There  are  three  different  reduction  schemes  included  in  the  algorithm:  non- 
K-vertex  deletion,  mdegree-one  K-vertex  contraction  and  series  and  parallel 
reductions.   However,  the  last  three  of  these  interact  directly. 

The  non-K-vertex  deletion  steps  (1)  and  (2)  in  REDUCT.  requires  0(|£|) 
operations  in  the  worst  case.  For  example,  consider  a  complete  acyclic  graph 
which  has  only  one  K-vertex,  the  source.  Here,  one  vertex  at  a  time  would  be 
deleted  from  the  network,  but,  in  the  process,  every  edge  would  be  examined 
exactly  once  until  a  single  isolated  vertex  remained.  The  initial  queue -building  is  cf 
no  consequence  since  it  is  0(  |  V\). 

Steps  (3)  through  (5)  constitute  the  other  three  reductions.  Queue-building 
again  may  be  disregarded  since  it  is  at  an  0(1  V\)  operation.  Initial  parallel 
reductions  may  require  0(|j?|)  operations  using  a  bucket  sort  technique  to  check 
for  edges  with  common  end  vertices.  The  central  step  of  the  algorithm,  step  (5). 
requires  0(|7|2)  time.    This  is  true  since  (1)  at  most   \V\   series  reductions  or 

30 


contractions  can  ever  occur  (2)  each  such  reduction  requires  one  0( ;  V\)  checK  for 
a  newly  created  parallel  edge,  and  (3)  once  identified,  actually  carrying  out  the 
parallel  reduction  requires  only  constant  time.  Overall  then.  REDUCT  requires 
0(|  V\2)  time. 

There  are  two  different  decomposition  schemes  used,  cut-vertex 
decomposition  and  neck-vertex  decomposition.  The  complexity  of  the  cut-vertex 
algorithm  is  0(|i?|)  since  it  is  based  on  the  0{\E\)  depth-first  search.  The 
complexity  of  the  neck-vertex  algorithm  Will  be  0(  |  V-K\  \E\)  since  finding  neck- 
vertices  is  based  on  repeated  0(  |  E  |)  searches  on  D-v  for  each  non-K-vertex  v .  In 
the  worst  case  then,  this  decomposition  requires  0(  |  V\  \E\)  time.  However,  if  we 
fix  the  number  of  non-K-vertices  as  described  in  chapter  4,  then  the  complexity  of 
neck-vertex  algorithm  will  effectively  be  0(  \E  |). 

It  follows  from  the  above  discussion  that  the  total  complexity  of  the  reduction 
and  decomposition  algorithms  is  0(!7|l£'|)  normally,  and  Od^l2)  when  non-K- 
vertices  are  fixed. 

After  reductions,  computing  SKT  reliability  requires  0(  \E\2'[V~K\)  time.  This  is 
obviously  true  since  Z^V~K^  combinations  of  non-K-vertices  must  be  produced  and 
an  0(j  £  |)  SAT  computation  made  for  each  combination.  "tVhen  the  number  of  non- 
K-vertices  is  fixed,  the  computation  becomes  0(|£*|).  The  entire  algorithm 
therefore  requires  0( \E\2^V~K'^)  time  in  general  but  only  0(  |  V\z)  time  with  non-K- 
vertices  fixed. 
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VI.      COMPUTATIONAL  EXPERIENCE 

In  order  to  test  the  efficiency  of  the  SKT  reliability  algorithm,  we  used  two 
kinds  of  networks,  "complete"  acyclic  networks,  and  acyclic  "street"  networks.  A 
complete  acyclic  network  is  a  graph  D{V,E)  such  that  V=\vl,V2,  •  •  ■  ,vn\  and 
e=(vitVj)eE  if  and  only  if  i<j.    Any  such  network  with  n  vertices  has  exactly 

— »—■ — *- edges.   A  network  is  called  anmxn  "street"  network  if  (1)  its  vertices  can 

be  laid  out  in  an  m  rows  by  n  columns  rectangular  grid,  (2)  for  all  vertices  except 
those  in  the  last  column  there  is  an  edge  incident  out  each  vertex  directed  to  the 
right,  and  (3)  for  all  vertices  except  those  in  the  first  row,  there  is  an  edge  incident 
out  of  each  vertex  directed  upward.  Whenn=7n  case  the  street  network  is  called 
"square"  (See  Figure  6.1). 

In  the  complete  network  we  assigned  vl  to  be  the  source  and  fixed  vn  as  a  PC- 
vertex.  This  insures  that  all  edges  and  vertices  are  relevant.  Other  K-vertices  were 
chosen  randomly  as  required.  In  the  street  network,  all  corner  vertices  were  fixod 
as  K-vertices,  the  lower  left-hand  corner  vertex  assigned  as  the  source  and  ail 
other  K-vertices  chosen  randomly. 

A      SOURCE-TO-ALL-TERMINAL  RELIABILITY  COMPUTATION 

As  described  in  chapter  4,  the  complexity  of  computing  SAT  reliability  is 
0(\E\).  Since  repeated  computation  of  SAT  reliability  is  the  key  to  computing 
general  SKT  reliability,  it  is  interesting  to  see  how  much  CPU  time  is  required  for  a 
single  SAT  computation.  Figure  6.2  shows  the  CPU  time  used  as  a  function  of  the 
number  of  vertices  in  the  two  types  of  test  networks.  |  V\  is  increased  from  15  to 
80  for  the  complete  network  and  from  9  to  144  for  the  square  street  network.  Of 
course,  as  a  function  of  \V\,   \E\   is  quadratic  for  the  complete  network  ranging 
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be  laid  out  in  an  m  rows  by  n  columns  rectangular  grid,  (2)  for  all  vertices  except 
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right,  and  (3)  for  all  vertices  except  those  in  the  first  row,  there  is  an  edge  incident 
out  of  each  vertex  directed  upward.  When  n=m  case  the  street  network  is  called 
"square"  (See  Figure  6.1). 

In  the  complete  network  we  assigned  vx  to  be  the  source  and  fixed  vn  as  a  K- 
vertex.  This  insures  that  all  edges  and  vertices  are  relevant.  Other  X-vertices  were 
chosen  randomly  as  required.  In  the  street  network,  all  corner  vertices  were  fixed 
as  K-vertices,  the  lower  left-hand  corner  vertex  assigned  as  the  source  and  all 
other  K-vertices  chosen  randomly. 

A      SOURCE-TO-ALL-TERMINAL  RELIABILITY  COMPUTATION 

As  described  in  chapter  4,  the  complexity  of  computing  SAT  reliability  is 
0(\E\).  Since  repeated  computation  of  SAT  reliability  is  the  key  to  computing 
general  SKT  reliability,  it  is  interesting  to  see  how  much  CPU  time  is  required  for  a 
single  SAT  computation.  Figure  6.2  shows  the  CPU  time  used  as  a  function  of  the 
number  of  vertices  in  the  two  types  of  test  networks.  |  V\  is  increased  from  15  to 
80  for  the  complete  network  and  from  9  to  144  for  the  square  street  network.  Of 
course,  as  a  function  of  [ Vj.   \E\  is  quadratic  for  the  complete  network  ranging 
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from  105  to  3160.  \E\  increases  linearly  as  a  function  of  \V\  for  the  street 
network  ranging  from  12  to  284.  From  the  figure,  it  is  clear  that  SAT  reliability  can 
be  computed  very  rapidly. 

B.  SOURCE-TO-K-TERMINAL  COMPUTATION 

Here  we  investigate  the  work  required  by  the  algorithm  for  computing  general 
SKT  reliability  and  the  efficiency  gained  by  using  reductions  and  decomposition. 
Initial  reductions  and  decomposition  require  0(  |  V\  \E  |)  time.  After  reductions  and 
decomposition,  the  complexity  of  computing  SKT  reliability  is 
0(max£  1 2   *  ^   |  |  E^  j  |)  where  C  is  the  set  of  separable  components  of  D  and  where 

icC 

Vi-Ki  is  the  set  of  non-K-vertices  in  component  i.  Consequently,  we  expect  an 
exponential  growth  in  the  worst  case  for  SKT  computations.  Figure  6.3  shows  the 
growth  in  CPU  time  as  the  number  of  non-K-vertices  is  increased,  with  and  without 
reductions  and  decomposition 

For  the  complete  network,  the  number  of  non-K-vertices  ranged  from  3  to  14 
out  of  a  total  20  vertices.  The  reductions  and  decomposition  do  not  reduce 
computation  much  here  since  the  complete  network  is  so  densely  connected.  On 
the  other  hand,  computation  times  are  significantly  reduced  when  reductions  and 
decomposition  are  used  in  the  street  network,  primarily  because  many  non-K- 
vertices  are  changed  to  K-vertices.  Notice  that  the  computation  time  required  by 
the  street  network  without  reductions  and  decomposition  is  almost  the  same  as  for 
the  complete  network.  This  shows  how  complexity  is  dependent  mostly  on  the 
number  of  non-K-vertices. 

C.  COMPUTATION  WITH  NON-K-VERTICES  FIXED 

If  the  number  of  non-K-vertices  is  fixed,  i.e.  \V-K\=c  for  some  constant  c, 
then  the  complexity  of  computing  SKT  reliability  is  0(  |  7':2)+0(c  \E\)  -  0(2C  \E\) 
which  is  O^Vl2).    Thus,  the  complexity  is  polynomial  instead  of  exponential.    Of 
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course,  if  reductions  and  decomposition  are  not  used,  the  complexity  becomes 
0{\E\)  which  is  better  than  0(]7j2).  However,  in  practice,  the  reductions  and 
decomposition  are  very  efficient  and  do  not  seem  to  increase  computation  time. 
Figure  6.4  shows  the  increase  in  computation  for  both  test  networks  as  \V\  is 
increased. 

D.      LOWER  BOUND  COMPUTATION 

As  described  in  chapters  4  and  5,  a  lower  bound  on  SKT  reliability  is  obtained 
when  the  normal  SKT  reliability  computation  is  truncated.  Two  heuristic 
arguments  were  given  in  chapter  4  indicating  that  this  lower  bound  should  be  fairly 
good.  Figure  6.5  shows,  for  actual  computations,  how  the  lower  bound  approaches 
the  exact  network  reliability  very  quickly.  The  test  networks  were  assigned 
different  values  of  pt  to  show  that  the  accuracy  of  the  bound  is  not  much  affected 
by  Pi-  It  appears  that  a  CPU  time  limit  can  be  used  successfully  for  truncating  SKT 
computation  to  obtain  a  good  lower  bound  when  exact  computation  is  not  feasible. 
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VII.     CONCLUSION 


This  thesis  has  explored  the  source-to-K-terminal  reliability  analysis  problem 
in  acyclic  directed  networks.  Vfe  have  shown  that  the  polynomial-time  reductions 
and  decompositions  reduce  an  exponentially  complex  problem  to  a  polynomially 
complex  problem  under  certain  conditions.  Even  when  these  conditions  do  not 
hold,  surprisingly  large  networks  can  be  analyzed.  The  theoretical  and 
computational  results  obtained  here  can  be  used  by  researchers  in  analyzing  the 
reliability  of  directed  acyclic  networks  and,  in  some  cases,  for  the  partial 
reliability  analysis  of  general  directed  networks. 

We  assumed  that  all  vertices  work  perfectly  and  that  edge  reliabilities  are 
constant  while  the  network  system  is  alive.  These  assumptions  may  not  be 
acceptable  in  a  real-world  situation.  It  may  be  necessary  to  treat  vertex 
reliabilities  explicitly  and,  in  some  way,  handle  component  reliabilities  changing  as 
a  function  of  time.  This  will  complicate  reliability  analysis  but  does  not  make 
analysis  infeasible. 

The  extension  to  unreliable  vertices  is  not  hard.  All  K-vertices  must  T.vork  if 
the  system  is  to  work  and  thus,  for  unreliable  K-vertices,  we  obtain  network 
reliability  as  usual  and  then  multiply  by  the  probability  that  all  K-vertices  function. 
Unreliable  non-K-vertices  may  be  handled  in  several  ways.  One  of  the  simplest 
ways  is  to  split  each  unreliable  non-K-vertex  into  two  perfectly  reliable  vertices 
with  an  unreliable  edge  between  the  two  vertices  with  a  failure  probability  equal  to 
the  failure  probability  of  the  original  vertex.  All  edges  incident  into  the  original 
vertex  are  made  incident  into  the  tail  vertex  of  the  new  edge  and  all  edges  incident 
out  of  the  original  vertex  are  made  incident  out  of  the  new  head  vertex.    However, 
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it  is  also  possible  to  redefine  the  reduction  and  decomposition  schemes  and 
redefine  the  general  reliability  expression  to  directly  handle  unreliable  non-K- 
vertices.  For  example,  the  series  reduction  can  include  the  vertex  reliability  so 
that  pc  =  paPbP(vi)-  Also,  the  term  in  the  general  reliability  expression  associated 
with  "no  edges  in  E^.,m  work"  can  be  replaced  with  a  term  which  expresses  "no 
edges  into  Vi  work  or  v.:  fails,  and  no  edges  into  Vj  work  or  Vj  fails,  and  ...  except 
edges  going  between  pairs  of  vertices  in  \v^,Vj,  •  ■  ■  ,vk\.  This  is  not  hard  to 
compute. 

Reliability  changing  as  a  function  of  time  can  be  handled  in  one  of  two  ways. 
The  easiest  way  is  to  just  input  component  reliabilities  at  different  times  and 
compute  network  reliability  at  these  times.  Interpolation  can  be  used  to 
approximate  values  between  the  times  for  which  reliability  is  actually  evaluated. 
Another  possibility  is  to  produce  a  symbolic  reliability  expression  that  can  be 
repeatedly  evaluated  by  a  computer  at  different  times,  thus  avoiding  repeated 
reductions,  decompositions,  etc.  This  might  require  large  amounts  of  storage  but, 
if  reliability  can  be  calculated  at  all  in  a  reasonable  amount  of  time,  it  should  be 
feasible.  Once  the  price  is  paid,  computing  reliability  for  any  given  point  in  time 
should  be  very  efficient. 

As  modern  computer,  power  and  other  network  systems  become  larger  and 
more  complicated,  reliability  becomes  a  more  important  part  of  design  criteria. 
Unfortunately,  it  is  obvious  that  as  networks  become  more  complicated  they  also 
become  more  difficult  to  analyze.  Further  research  should  concentrate  on  ways  to 
streamline  computations  and  to  find  accurate  but  easily  computable  bounds  The 
method  described  in  chapter  4  for  computing  SKT  reliability  could  be  improved  if 
only  non-zero  terms  were  generated;  An  enumeration  scheme  similar  to  that  of 
Satyanaryana  and  Prabhakar  [1978]  should  help  in  this  respect.  The  lower  bound 
described  in  this  thesis  is  a  good  start  on  finding  accurate  and  easily  'computable 


36 


bounds,  but  every  lower  bound  should  be  accompanied  with  an  upper  bound  so  that 
accuracy  can  be  checked.  Finding  a  good  upper  bound  should  be  a  priority  for 
future  research. 
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APPENDIX  A 


<a) 
v  a  cut-vertex 


(b) 
v  not  a  cut-vertex 


Neck  Vertex  to  K-vertex  Reduction 
<Figure  3.1> 


Decomposition 
<Figure  3.2> 
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Moving  Edges  to  Source 
<FUure  3.3> 


All  Terminal  Reliability  Rsv{D\)  =  RsvWs) 
<Figure  4.1> 
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Completely  Reducible  Network  by  Reduction  Schemes 

<Fiffure  5.1  > 
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Multilist  Data  Structure 
<?igure  5.2> 
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4x4  Square  "Street"  Network 
<Figure  6.1> 
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Source-To-All  Terminal  Network: 
<  Figure  6.2  > 
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ELAP5ED  CPU  TIC  WITH  AM)  WTTHOUT  REDUCTION.  DECOMPOSmON 
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Souce-To-A'  Terminal  Network 
<  Figure  6.3  > 
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COMPUTER  CPU  TIC  WITH  FDSD  NON-K-VERTKXS  10 
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So'drce-To-K  Terminal  with  Fixed  Non-ZT- Vertices 
<  Figure  6.4  > 
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RQUBIUTY  INCREASE  8Y  TIME  IN  TOTAL  20  VERTICES 
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Network  Reliability  Increment 
<  Figure  6.5  > 
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